Draft version February 5, 2008 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



O 
O 



SUPERNOVA REMNANTS IN THE MAGELLANIC CLOUDS. VII. INFRARED EMISSION FROM SNRS 

R. M. Williams, Y.-H. Chu, and R. Gruendl 

University of Illinois at Urbana-Champaign, 1002 W. Green St., Urbana, IL 61801 USA 
Draft version February 5, 2008 

ABSTRACT 

We have used the instruments on the Spitzer Space Telescope to study the Large Magellanic Cloud 
supernova remnants (SNRs) N11L, N44, N49, N206, N63A, and N157B. The two large SNRs N44 and 
N206 were not detected in any IRAC or MIPS wavebands; the remainder were detected at one or 
more wavelengths. In particular, the SNRs N49 and N63A each had features that were evident in all 
available IRAC and MIPS bands. Each of these two also displayed faint limb emission in the MIPS 
24 /im band only. IRS spectra obtained for the N49 SNR showed a number of prominent lines, with 
little continuum contribution. We therefore suggest that N49, and possibly N63A, are dominated by 
line emission, with thermal emission from hot dust being at most a secondary component. 
Subject headings: ISM: supernova remnants - ISM: individual (N11L, N44, N49, N206, N63A, N157B) 
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1. INTRODUCTION 

Supernova remnants (SNRs) are expected to produce 
infrared (IR) emission from such mechanisms as atomic 
and molecular line emission (dominated by ground-state 
fine structure lines), free- free emission from hot gas, syn- 
chrotron emission, and thermal continuum emission from 
dust heated by collisions with the post-shock gas. Of 
these, we expect onl y negligible contribut ions from free- 
free emission fe.g. . iFis chera et al.ll2002fl or from syn- 
chrotron emission |Grah^ni^t^iiTll^87a> . except where 
a pulsar-wind nebula (PWN) may be present. Thus, the 
IR emission is thought to be dominated by line emission 
and/or thermal emission from dust. The dust emission 
is of particular interest to the study of SNRs, as super- 
novae are thought to produce a s ignificant fraction of 
dust in the ISM (e.g., iDweld 119981 Dust emission, in 
turn, is an important factor in cooling the hot plasma 
within S NRs through inelastic collisions with elec trons 
and ions llOstriker fc Silldll973l ISmith fc Coxll200l . 

Infrared emission was detected with the Infrared As- 
tronomy Satellite (IRAS) from roughly a third of the 
known Galactic SNRs, demonstrating the difficulty in de- 
tecting many of these objects in the confusion with other 
IR sources in the ISM. Because of this confusion, most 
of the SNRs to be observed were comparatively bright, 
nearby objects l)ArendtJll989|) . Similarly, IRA S observa- 
tions of the Large Magellanic Cloud (LMC), including 
pointed observations toward nine known re mnants, only 
detected four SNRs ijGraham et al J 119874 - Using the 
more rece nt G LIMPSE surve y from the Spitzer Space 
Telescope, iReach et al.l ()2006(l clearly detected 18 of 95 
Galactic SNRs known to be within the fields. 

Spitzer observations are sensitive to many forms of in- 
frared emission from SNRs, depending on the instrument 
and wavelength b and. The Infrared Spectrograph (IRS; 
iHouck et alJl2004)) is ideal for the st udy of line emissi on. 
The Infrared Array Camera fIRAC: lFazio et alJl2004|) is 
expected to include emission from polycyclic aromatic 
hydrocarbons (PAHs) and very small grains (VSGs); H 2 
lines; and a number of atomic lines including Bra 4.1 
/im, [Fe II] 5.3 /rnr and [Ar II] 7.0 /im. The Multi- 



band Imaging Photometer for Spitzer fMIPS: lRieke et all 
2004) bands are also expected to include some atomic line 
emission, including [Fe II] 24.5 /im, [O I] 63.2 /im, and 
[C II] 157 /im; additionally, we expect significant thermal 
emission from collisionally heated dust. 

Dust emission from SNRs is of particular interest. It 
has been suggested that supernovae (SNe) a re one of the 
principal sources of dus t production (e.g., iDwekl Il998t 
ITodini fc Ferraral l200ll iDunne et all 12003) . Theoreti- 
cally, models have predicted as much as 0.2-4 M of 
dust production in a typical Type II sup ernova (e.g., 
ITodini fc Ferraral 120011: IDunne et all 12003) . Certainly 
spectroscopic and photometric observations have sug- 
gested dust fo rmation in at le ast five supernovae, includ- 
ing SN1987A l|Datizigerll2005T) . 

Finding this freshly produced dust in SNRs, however, 
has proven problematic. Infrared fluxes attributed to 
warm (~ 100-200 K) dust have been observed in youn 
remnants such as Cas A Tycho, and Kepler (e.g.. lBra~ 
H987tlArendt et alJll99a iHines et alJl2004|) . but the in- 
ferred dust mass is generally several orders of magnitude 
less than that predicted. Furthermore, it is often unclear 
whether this dust was largely produced in the SN or has 
since been swept up from the surrounding medium. 

Dust at cooler temperatures (< 25 K) has also been 
detec ted toward rem nants such as Cas A lIDunne et alJ 
12003) and Kepler ijMorgan et all 12003) using sub- 
millimeter continuum observations. Substantially 
greater dust masses (> 1 M Q ) have been deduced from 
these observations, but again, how much of this dust is 
ejecta from the recent SN is unknown. Worse, not all 
of the observed dust may be p hysically associated with 
the SNRs themselves. In fact, iWilson fe Batrlal l|2005f) 
suggest that "at least one-half of the "cold" dust emis- 
sion from the direction of Cas A is actually foreground 



emission. 



2. OBSERVATIONS 

We obtained proprietary Spitzer observations for re- 
gions which include the LMC SNRs N63A, N11L and 
SNR 0523-67.9 in N44, as well as archival Spitzer data 



. for SNRs N49, ,SNR 0532-71.0 in N206, and N157B. The 

Electronic address: rosanina@astro.uiuc.edu. chu@astro.uiuc.edu, gruendl@astro.uiuc.edu , . „ , , m a n j ±. 
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TABLE 1 
Spitzer Datasets used 



SNR 


Inst. 


Mode 


Prgrm 


AOR 


PI 


Date 


N11L 


IRAC 


5 x30s Man 


3565 


11 1 71 840 

111 / lu^u 


Chu 


2004 Nov 30 


N11L 


MIPS 


Sccin rVleip 


oovo 


111 77798 
111 / t I zo 


Chu 


900^ A/Tar 7 


N44-SNR 


IRAC 


5 x30s Map 


3565 


11172352 


Chu 


2005 Mar. 28 


N44-SNR 


MIPS 


Scan Map 


3565 


11177984 


Chu 


2005 Apr. 7 


N49 


IRAC 


10 x2s Map 


124 


8152064 


Gehrz 


2004 May 5 


N49 


MIPS 


24/tm Phot 


124 


8151808 


Gehrz 


2004 May 26 


N49 


MIPS 


70ftm Phot 


124 


8791040 


Gehrz 


2005 May 26 


N49 


IRS 


Stare 


124 


6586112 


Gehrz 


2004 May 31 


N206-SNR 


IRAC 


3 x12s Map 


1061 


6063104 


Gorjian 


2003 Nov. 21 


N206-SNR 


MIPS 


Scan Map 


717 


7864320 


Riekc 


2003 Nov. 24 


N63A 


IRAC 


5 x30s Map 


3565 


11173888 


Chu 


2004 Dec. 16 


N63A 


MIPS 


Scan Map 


3565 


11178496 


Chu 


2005 Mar. 8 


N157B 


IRAC 


3 x12s Map 


63 


4379904 


Houck 


2004 Jan. 12 


N157B 


IRAC 


3 x12s Map 


1032 


6056960 


Branch 


2003 Nov. 06 



were processed by the Spitzer Science Center (SSC) to 
produce Basic Calibrated Data (BCD) files. 

Images with IRAC were obtained for all six SNRs. The 
observations were carried out in mapping mode for the 
3.6, 4.5, 5.8, and 8.0 /zm bands. All observations ex- 
cept that for N49 were made with high dynamic range 
and cycling dither; the N49 images were taken with a 
5-Gaussian dither. The exposure times for each map are 
included in Table 2] We used post-BCD mosaic images 
resulting from the standard SSC reduction pipeline (ver- 
sion Sll). 

MIPS scan maps were available for SNRs N63A, N11L, 
N44 (SNR), and N206 (SNR) in the 24, 70, and 160 /im 
bands. All of the observations were taken at a Medium 
scan rate, with a scan leg length of 0?5 and a scan step 
of 148". Spatial coverages were 50'x20' for N11L, N44 
(SNR), and N206 (SNR), and 50'x 15' for N63A. For N49, 
observations in the MIPS Photometry mode were avail- 
able in the 24 and 70 /jm bands. The 24 /im observation 
took four 3 s observations of a small field, while the 70 
/im observation took one 10 s observation of a large field. 
We used the standard post-BCD images for our study. 

N49 was observed with the IRS in staring mode. Low- 
resolution spectra were obtained in the 5.2-8.7 /im (SL2) 
and 7.4-14.5 /im (SL1) bands; high-resolution spectra 
were obtained in the 9.9-19.6 /im (SH) and 18.7-37.2 /im 
(LH) bands. For each instrument, a spectrum is taken at 
two "nod" positions, intended for on- and off-source com- 
parisons. In this case, both positions fall on similar fila- 
mentary regions in N49. IRS spectra were extracted and 
calibrated using the SSC SPICE package. The flux cali- 
bration for these spectra is based on SSC established con- 
version factors from electrons s _1 to Jy, corrected by fac- 
tors determined from the spectra of standard stars. No 
background subtraction was performed, as both "nod" 
positions lie along emitting material within N49. 

3. RESULTS 
3.1. Imaging 

There is a great deal of variation as to which SNRs 
can be seen in each of the IRAC and MIPS wavebands 
(above the level of the nearby H II complexes). Two of 
these SNRs, N44 and N206, are not detected in any band; 
two others, N11L and N157B, are faintly visible at 4.5 
/im. N63A and N49 both have interior concentrations of 



material that are evident at all available wavelengths; in 
both cases, however, the SNR limb is distinct at 24 /im 
but undetectable in other wavebands. 

Flux densities for these objects in each waveband are 
listed in Table [21 When the object is not detected in a 
given waveband, an upper limit for a detection 2a over 
the background is given instead. These fluxes and upper 
limits are preliminary figures and should be treated with 
caution. For one, the SNRs are generally superposed on 
diffuse emission from the nearby H II complexes. To ac- 
count for this and other background emission, we found 
the median level of emission immediately surrounding 
each SNR, and subtracted this from the measured flux 
values. The uncertainties in the background estimates 
suggest that the listed fluxes are uncertain by 20-30%. 
In addition, there are uncertainties in the photometric 
calibrations of both IRAC and MIPS, ranging from at 
or below 10% for IRAC and for MIPS 24 /mi, to 20% 
for MIPS 70 and 160 /im. Because these objects are ex- 
tended, the IRAC fluxes are multiplied by the correction 
factor for an "infinite" aperture as recommended in the 
IRAC Data Handbook 1 : 0.94, 0.94, 0.63 and 0.69 for 3.6, 
4.5, 5.8 and 8 /im, respectively. The SNRs are discussed 
individually below. 



3.1.1. N11L 

N11L (SNR 0454-68.5) is on the periphery of the Nil 
H II complex, framed by a 150 pc filament at the edge of 
that complex. The SNR itself shows a clear shell struc- 
ture in optical emission lines, with filaments projecting 
beyond the boundary of t hat shell. Its optical e xpansion 
velocity is ~ 200 km s" 1 ijWilliams et al.lll999|) . 

N11L was close to the edge of the region scanned in the 
Nil observation; it was only partially covered in the 3.6 
and 5.8 /im IRAC bands, and was at the extreme edge of 
the three MIPS bands. No emission from the SNR was 
distinguished in any of these bands. Likewise, the 8 /im 
IRAC band showed no emission above the background 
level for this SNR. In the 4.5 /im IRAC band, however, 
the SNR shell is faint but distinct, with a similar mor- 
phology to that seen in Ha (Fig. [TJ. 

1 http: / /ssc. spitzer. caltech.edu/irac/dh/iracdatahandbook2. 0.pdf, 
Table 5.7 
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TABLE 2 
Flux Density Estimates 



Band 


N11L 


N44 


N49 


N49-limb 


N206 


N63A 


N63A-limb 


N157B 


H 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mjy) 


(mJy) 


(mjy) 


3.6 




< 74 


32 


<8.8 


< 96 


37 


< 6.8 


62 


4.5 


1.5 


< 51 


41 


<8.2 


<70 


32 


< 5.9 


87 


5.8 




< 86 


130 


<26 


< 140 


130 


< 17 


< 140 


8.0 


< 8.5 


< 190 


180 


<30 


< 250 


300 


< 45 


< 210 


24 




< 744 


1500 


66 


< 300 


2300 


760 




70 




< 5100 


10,200 


< 433 


< 3200 


7700 


< 2100 




160 




< 12,500 






< 2200 


5800 


< 2300 




Area" 


0.73 


8.80 


1.11 


0.40 


8.31 


1.31 


0.92 


1.20 



Note. — "<" indicates that the SNR is not detected in this waveband. The number 
given is an upper limit, giving the flux required for a detection 2c above the background. 

a Area on the sky (of SNR or portion thereof) for which flux densities are estimated, given 
in square arcminutes. The N49-limb area includes only the eastern limb, without any of the 
bright filaments; the N63A-limb area excludes the three-lobed nebular emission; and the 
N157B area excludes the southern dust clouds. 




Fig. 1.— Images of N11L with IRAC (a) 4.5 and (b) 8.0 ^m; and 
with the Cerro Tololo Inter-American Observatory (CTIO) Curtis 
Schmidt Telescope in (c) Ha and (d) [S II]. 




Fig. 2.— Images of N44 with IRAC (a) 5.8 and (b) 8.0 A*m; 
MIPS (c) 24 ^im; and the CTIO Curtis Schmidt Telescope in (d) 
Ha. 



3.1.2. SNR 0523-67.9 in NU 

The N44 SNR is to the northeast of the N44 H II re- 
gion, overlapping with the H II region to the west and 
south ijChu et al.i lT9931. The SNR has a bi-lobed struc- 
ture, possibly due to an intersecting dust lane. It shows 
an i rregular expansio n, with velocities of up to 150 km 
s" 1 <|Kim et al.111998^ . 

The SNR in N44 is, in general, conspicuous by its ab- 
sence in the IRAC and MIPS images. There seems to 
be no infrared emission at Spitzer wavelengths following 
the SNR shell seen in the optical and X-ray regimes. A 
large cloud seen in all bands overlaps the SNR south- 
western limb; an "arm" of that cloud seen clearly at 5.8, 
8.0 and 24 /im protrudes across the SNR between its two 
lobes. This "arm" , seen at the IR wavelengths which are 
most sensitive to dust emission, is coincident with the 
dust lane causing obscuration in optical images. Some 
5.8 and 8.0 /jm emission surrounds the shell boundary, 



but aside from the intersecting "arm" , the emission in 
these bands is notably lower over most of the optical ex- 
tent of the SNR (Fig.|2J). While it is quite possible that 
some of the dust "arm" is shock-heated by the SNR be- 
hind it, the "arm" appears to be part of a much larger 
IR emission feature that extends well south of the SNR 
on the sky. 

3.1.3. N49 

N49 (SNR 0525-66.1) is unusual for this sample, in 
that it is not situated in or near an H II region. How- 
ever, it does lie nea r a molecular clou d, with which it 
may be interacting l)Banas et alJ 119971). X-rav and ra- 
dio maps show a com plete SNR shell l|Park et alJl2003t 
Dick el fc Milndll998j) . while optical maps largely high- 
light a co mplex of filaments c overing the eastern side of 
the SNR llBilikova et alJ|2006D. It was noted as an IRAS 
source bv lGraham et all <|1987a|) . 




4.fi Ulli 



Li 



BiiIjL btmibii IMWI 

Fig. 3.— Images of N49 with IRAC (a) 3.6 (b) 4.5 (c) 5.8 (d) 
8.0 fim and MIPS (e) 24.0 and (f) 70 fim. The E and W in panel 
(a) label the East and West sides of the SNR. 



Unlike most of the remnants discussed here, N49 shows 
clearly in all four of the IRAC bands and the 24 and 70 
/zm MIPS bands (Fig. 0. Observations of this SNR in 
the 160 /jm band are not yet available in the archive. 
Regions that are bright in the IRAC bands correspond 
well to brighter regio ns of emission at optical wavelengths 
(Bilikova~elall2006), even showing some of the same fil- 
amentary structure (Fig. @J- The 70 /im image shows 
very similar structure, albeit at lower resolution. At 24 
/im, however, in addition to these bright emission re- 
gions, one can also detect faint emission from the rest of 
the SNR as seen in the X-ray and radio regimes, with a 
morphology very similar to that in radio. In particular, 
the western limb of the SNR, which has no detectable 
optical emission, is seen faintly but distinctly at 24 /jm. 
In Table [3 we estimate the flux density from the western 
limb only, as well as that from the entire SNR. 

3.1.4. N63A 

N6 3A (SNR 05 35-66.0) lies within the N63 H II re- 
gion It was noted as an IRAS source by 
iGraham et alJ l)1987a|) . While radio and X-ray obser- 
vations of the SNR show a complete shell l)Dickel et al.l 
I19931 1 Warren et al.ll2003T) . most of this shell is not seen 
in optical emission lines. The bright optical emission is 
confined to a three-lobed st r uctur e on the western side 
of the SNR. ILevenson et alJ l|1995f) showed that the two 
eastern lobes of this nebula have high [S II] /Ha ratios, 
indicative of shock ionization, while the optical spectrum 
of the western lobe is more consistent with the photoion- 
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Fi g. 4. — Images o f N49 with (a) the Chandra X-ray Observa- 
tory IPark et al.l200 3T): (b) the A ustralia Telescope Compact Array 
at 13 cm iDickcl fc Milnelll99SI) : and the Hubbl e Space Telescop e 
(HST) WFPC2 in (c) Ha and (d) [O III] IBilikova et aT] 120061) 
The c ontours in panel (a) are taken from Fig. 5 of TBanasi'ltHIrt] 
119971) and show a local CO peak; these contours are at 1.0, 1.3, 
1.7, 2.0, 2.3 and 2.7 K km s" 1 . 



TABLE 3 
Flux Densities for N63A 



'Lobes" 



Band 


N shock 


S shock 


W phot 


(pm) 


(mJy) 


(mJy) 


(mJy) 


3.6 


4.6 


11 


11 


4.5 


4.7 


9.7 


9.0 


5.8 


13 


35 


30 


8.0 


30 


93 


87 


24 


280 


370 


370 


Area 


0.044 


0.054 


0.036 


I(3.6)/I(5.8) 


0.35 


0.31 


0.37 


I(4.5)/I(8) 


0.16 


0.10 


0.10 



Note. — The three lobes are referred to as 
"N shock" for the northeastern shocked lobe, "S 
shock" for the southeastern shocked lobe, and 
"W phot" for the western photoionized lobe. Ar- 
eas are in square arcminutes. 



ization in H II regions. 

The optically-bright nebula is also bright in the IRAC 
and MIPS wavebands (Fig.0. In the IRAC and 24 fj,m 
MIPS bands, where the resolution is sufficient to dis- 
tinguish between the lobes, the western lobe is notably 
brighter than the two eastern lobes. In the IRAC bands, 
this three-lobed region is the only feature of N63A to 
appear clearly above the background emission. (A few 
bright optical knots elsewhere in the SNR may have some 
very faint corresponding emission in the IRAC images, 
but that emission is very close to background levels.) 
This also appears to be the case for the 70 and 160 ^m 
observations, although the lower resolution of the latter 
images makes it more difficult to distinguish emission 
associated with the three-lobed nebula from the rest of 
the SNR. The flux densities for these lobes are given in 
Table □ 
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Fig. 5.— Images of N63A with (a) IRAC 8.0 /im; (b) MIPS 
24.0 Atm; (c) the H ST WFPC2 in Ha, and (d) smoothed X-ray 
data from Chandra I Warren ct al. 2003). A background galaxy is 
marked with a circle in (a). 



The 24 jjbm MIPS band, on the other hand, shows 
clear emission over the entire X-ray/radio shell of the 
SNR. The emission is evenly distributed over most of the 
face of the SNR, aside from the aforementioned three- 
lobed structure. A slightly brighter spot on the east- 
ern limb corresponds to an area with [S Ill-b right knots 
and enhanced X-ray emission (Fig. |SJ also IChu et all 
Il999h . The 24 /im image also shows faint loops to the 
southwest, projecting beyond the boundary of the SNR. 
These correspond we ll to soft X-ray "crescents" noted by 
I Warren et alJ l)2003j) . With the exception of the three- 
lobed nebula (seen in X-rays as a heavily absorbed re- 
gion), there is a strong correspondence between the 24 
/im emission and the X-ray emission seen by Chandra. 
The flux density of this 24 /im emission, as with N49, is 
given separately in Table as well as that of the SNR 
as a whole. 

3.1.5. SNR 0532-71.0 in N206 

The SNR lies just outside the eastern rim of the N206 
H II region, with some overlap toward the south of 
the SNR. It was obs erved by IRAS but not detected 
ijGraham et alJll987a|) . It has an estimated expansion 
velocity of 200 km s _1 . An elongated radio and X- 
ray feature within the remnant, associated with an X- 
ray point source, is th ought to be a pulsar-wind nebula 
( Willi am's et alJl2"005j) . Spitzer da ta for the neighborin g 
H II region have been analyzed bv lGoriian et alJ {2004), 
who note that the region is rich in PAH emission except 
for the northeast side, nearest the SNR. 

As with the SNR in N44, the SNR in N206 does not 
show emission easily distinguishable from that of the sur- 
rounding H II region (Fig. EJ. An extension from the 
H II region overlaps the southern limb of the SNR, but 
does not seem to be correlated with the emission of the 
SNR at other wavelengths. Spitzer images of the N206 
H II region, includin g the SNR, have been published by 
iGoriian et all l)2004|) . 




Fig. 6.— Images of N206 with IRAC (a) 5.8 and (b) 8.0 /im; 
MIPS (c) 24 /im; and the CTIO Curtis Schmidt Telescope in (d) 
Ha. 



3.1.6. N157B 

N157B (SNR 0538-69.1) adjoins the 30 Doradus giant 
H II region. It is one of the two LMC remnants for which 
the presenc e of an associate d pulsar has been directly 
confirmed. IChu et all l|1992|) used optical imaging and 
spectroscopy to determine a maximum SNR extent of 
about 90" x 70", with patchy optical emission over its 
face and no clear shell structure. Dark clouds appear to 
be obscuring the southern portions of the SNR. 

No MIPS data are yet available for N157B. In the 
IRAC bands (Fig. UJ , the brightest feature is a "knot" on 
the southern side of the SNR, with an embedded point 
source. This bright IR region corresponds well to a region 
of lower Ha emission. However, there is also fainter emis- 
sion at 3.6 and 4.5 /im north of this dust cloud. At both 
these wavelengths, the emission is somewhat enhanced 
(compared to the background) within the ellipse that en- 
compasses most of the bright optical emission (Fig. 0;). 
In particular, a pair of curved filaments to the east of the 
SNR are seen clearly, particularly at 4.5 /im. At 5.8 and 
8.0 /im, any emission from the SNR is confused with the 
extensive foreground/background emission in the region. 

The faintness of the infrared emission is particularly 
surprising when one considers the pulsar-wind nebula 
within this SNR. Radio observations have suggested a 
spectral index of a = -0.19 for the PWN [S v oc v a ; 
lLazendic et alJ feOOflj) . Extrapolating from this radio 
spectrum, we would expect synchrotron contributions 
from the PWN of 320-360 mJy in the 3.6-8 /tm range, 
particularly at the location of the PWN itself as seen in 
radio and X-rays. Instead, the infrared emission, where 
detectable above the background, has much lower esti- 
mated flux densities and seems to follow the more promi- 
nent optical filaments rather than the radio/X-ray elon- 
gated PWN structure. 

Of course, it is common for PWNe to exhibit a "bro- 
ken" power-law spectrum, wherein the spectral index 
steepens as one goes to shorter wavelengths. The spec- 
tral break is often attributed to "aging" by synchrotron 
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Fig. 7.— Images of N157B with IRAC (a) 3.6 (b) 4.5 (c) 5.8 
(d) 8.0 /an; and with the CTIO Curtis Schmidt Telescope in (c) 
H« and (f) [S II]. The optical boundary of the SNR is shown by a 
dotted ellipse in (e). 



losses of a population of accelerated particles initially 
distributed with a power-law spectrum. In N158A (SNR 
0540-69.3), the other LMC SNR to contain a PWN, 
this spectral break occurs at ultraviolet wavelengths; the 
7—15 /im fluxes measured from ISO CAM are consis- 
tent with an extrapolation from the radio spectral in- 
dex {Gallant fc Tuff si 11999). In N158A, also, the spa- 
tial distribution of the infrared emission with ISO closely 
matches that of the radio emission. Recent Spitzer obser- 
vations of N158A, however, suggest that its IRAC-band 
emission follows the extrapolation of a simple power- 
law from optical wavelengths rather than radio - that 
is, "above" the break in the broken power-law spectrum 
{Reynolds et al.ll2005ft . If the spectral break in N157B is 
similarly located at wavelengths longer than 8 //m, the 
resulting lower flux expected from synchrotron emission 
may be more consistent with our non-detection. 

3.2. IR Spectroscopy of N4-9 

We examined the available low- and high-resolution 
IRS spectra for the SNR N49 at both "nod" positions. 
The placement of the slits is shown in Fig. [5] Spectra 
for the first nod position are shown in Fig. |5J The most 
prominent lines are given in TableQ] with their measured 
central wavelengths and intensities. The line intensities 
are given in terms of both peak flux density ( Jy) and line 
flux (erg cm~ 2 s _1 ). Uncertainties in the flux calibration 
(for point sources) are estimated at la values of 20% for 
the low-resolution and 30% for the high-resolution obser- 




FlG. 8. — Positions of the IRS slits superposed on the Hubble 
Space Telescope (HST) WFPC2 Ha image of N49. From shortest 
to longest, these are the Short-High Resolution, Long-High Reso- 
lution, and Short-Low Resolution slits. 



vations. However, the gain along the slit has not yet been 
fully calibrated, so additional errors may be introduced 
to the flux calibrations for our extended sources. There 
appears to be a significant flux "jump" between similar 
lines observed by the SL and SH modules; in part this is 
probably due to the lack of background subtraction, but 
there is also a known mismatch between the detectors 2 . 
Tentative identifications are given for each line based on 
wavelength coincidence with ionic and molecular lines 
known to be promi nent over the 5.2-38 wavelength 
range in SNRs fe.g.. lQliva et al.lll999albh . 

4. DISCUSSION 

Of our sample of 6 known SNRs within our available 
fields, we had 2 strong and 2 marginal detections in at 
least one waveband, along with the detection of a fea- 
ture possibly associated with the N44 SNR. For compar- 
ison, 18 Galactic SNRs were strongly detected (score 1) 
and 17 possibly de tected in confused fields (score 2) by 
IReach et alJ ((2006) , out of a sample of 95 located within 
the fields of t he Spitzer GLIMP SE survey. Note, how- 
ever, that the IReach et alJ l)2006ft search did not include 
MIPS data, so features which appear primarily at 24 /im, 
such as the outer limbs of N63A and N49, would not be 
seen in their survey. 

4.1. Non- detections 

The non-detections of the SNRs in N206 and, possibly, 
N44 are relatively unsur prising. Bo th of these SNRs have 
radii >20 pc (Table 151 Ichu et al]ll993T: IWilliams et all 
12005ft . while the other four SNRs have radii of ~ 8 - 16 
pc. Thus, surface brightness becomes more of an issue for 
these larger remnants, and may well put them under the 
upper limits for detection in most wavebands (Table |2J). 
The lack of detectable emission from the N206 SNR is 
particularly notable in light of the absence of PAH emis- 
sion in that area of the N206 H II region, as noted by 
IGoriian et all pOolft . 

We expect the environments of these SNRs to also have 
a strong influence on their detections (or lack thereof). 

2 IRS Data Handbook, http:/ /ssc. spitzer. caltech.edu/irs/dh/irsDH10. pdf 
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TABLE 4 

Lines and Line Fluxes from IRS data for N49 



Nod Position 1 



Nod Position 2 



Line 



A(/im) 


F„ a 


F b 


A(/im) 


F^ a 


F b 


ID 


Comments 


Short- 


-Low Resolution 










5.34 


0.011 


11 


5.34 


0.013 


8.2 


[Fe II] 


2nd order 


7.03 


0.016 


11 


7.00 


0.027 


18 


r A TTl 

[Ar 11] 


2nd order 


9.72 


0.033 


16 


9.72 


0.047 


19 


H 2 


(0,0)S(3) trans 


12.32 


0.034 


14 


12.32 


0.041 


14 


H 2 


(0,0)S(2) trans 


12.87 


0.189 


46 


12.87 


0.193 


56 


[Ne II] 


1st order 


Short- 


-High Resolut 


ion 








10.54 


0.132 


3.8 


10.50 


0.071 


0.8 


[S IV] 


noise 


12.30 


0.180 


3.9 


12.30 


0.153 


5.3 


H2 


(0,0)S(2) trans 


1Z.OO 


1.00 


47 


12.83 






[IN C 11J 


15.57 


1.15 


43 


15.57 


0.567 


20 


[Ne III] 




17.05 


0.333 


11 


17.05 


0.266 


7.9 


H 2 


(0,0)S(1) trans 


17.94 


0.665 


24 


17.94 


0.297 


13 


[Fe II] 


18.73 


0.376 


14 


18.73 


0.183 


8.6 


[S III] 




Long— 


High Resolution 








22.95 


0.839 


18 


22.95 


0.549 


12 


[Fe III] 




24.55 


0.247 


19 


24.55 


0.865 


20 


[Fe II] 
[O IV] 




25.90 


2.71 


55 


25.90 


2.75 


54 




26.00 


5.83 


125 


26.00 


6.28 


115 


[Fe II] 




33.47 


2.37 


49 


33.47 


2.62 


44 


[S III] 




34.84 


12.9 


263 


34.84 


14.3 


217 


[Si II] 




35.36 


2.27 


52 


35.36 


2.67 


46 


[Fe II] 




Note 


— The 


"comments" field includes (a) whether the line is seen in 



first or second order for the Short- Low spectra; (b) which specific transition 
of molecular hydrogen causes a given line; and (c) whether the line is 
notably affected by noise in that region of the spectrum. 

a Peak flux density in Jy 

b Line total flux in 10 — 14 erg cm -2 s — 1 




10.0 11.0 12.0 



13.0 14.0 15.0 16.0 
Wavelength, 



17.0 18.0 19.0 




Wavelength, |jm 




22.5 



35.0 



Wavelength, urn 
Tune 




11.0 12.0 
Wavelength, urn 



Fig. 9. — Spectra of N49 with IRS: (a) Short-High Resolution, (b) Long-High Resolution, (c) Short-Low Resolution, 2nd-3rd order, and 
(d) Short-Low Resolution, 1st order. 



To show substantial dust emission would require the pres- 
ence of dense clouds interacting with the SNR. There 
is no evidence for the presence of such dust clouds in 
the immediate vicinity of the N206 SNR; in N44, the 
optically-obscuring "arm" suggests at least some dense 
concentrations, but not sufficient for shock-heated dust 
to raise the SNR's IR surface brightness over the back- 
ground. 

4.2. Detections in 1 or 2 bands 

The detection of weak emission from N11L at 4.5 /im 
and N157B at 3.6 and 4.5 /xm are likely to be due to 
line emission from these remnants. Neither SNR is dis- 
cernable in the dust-sensitive 8.0 /im band; the shorter- 
wavelength bands, in general, are dominated by ionic and 
(espe cially) molecular line emission (e.g., iReach et al.l 
2006). In N11L, the 4.5 /im morphology closely traces 
the SNR limb and filament structures seen in optical 
emission lines. In N157B the 3.6 and 4.5 /im emission 
are enhanced over the region where optical filaments are 
brightest, and in particular match the morphology of a 
pair of bright optical filaments. The relative prominence 
of the 4.5 /im band, in both cases, suggests that Bra 4.1 
/im may be a strong contributor to the IR emission. 

4.3. Detections in multiple bands 
4.3.1. Origins of infrared emission in N^9 

N49 is bright at all of the wavelengths observed to date 
by Spitzer. This may be a result, in part, from N49's in- 
teraction with a molecular cloud to its southeast. In this, 
it is similar to the Galactic SNR IC443, which is also in- 
teracting with a nearby molecular cl oud, and is bright at 
a wide range of i nfrare d wavelengths llMufson et al.ll986t 
IBrann Rr. StwZttWMl Who et alJl2001l) . The TRS ohser- 
vation of N49 confirms the presence of numerous strong 
ionic lines, as well as several lines probably originating 
from H 2 transitions. The SNR does not appear to be 
dominated by PAH/VSG emission, as we would expect 
PAH-dominated emission to be strongest in the 8 /im 
band, while in fact the emission there is somewhat closer 
to background levels than at 3.6 or 4.5 /im. The lack 
of PAH emission could be accounted for by, e.g., PAH 
destruction by far-ultraviolet (FUV) radiation from the 
shock precursor of N49, which is known to ha ve fairly 
rapid expansion velocities l)Bilikova et alJl2~006D . 

We then consider whether thermal emission from dust, 
as well as line emission, makes up a significant part of 
the IR emission from N49. The emission in the IRAC 
bands only traces the brightest optical filaments, which 
also show up strongly in the MIPS 24 and 70 /im bands. 
These bright filaments most probably are the source of 
the emission measured in the IRS spectra, which indi- 
cate a number of strong emission lines. The continuum 
levels of the post-BCD spectra appear to be very close 
to zero, strongly suggesting a lack of substantial contin- 
uum emission from hot dust. We therefore argue that 
line emission is the dominant component of the mid/far 
infrared emission from the bright filaments which are the 
dominant source of emission from N49. While we cannot 
rule out additional contributions from hot dust, this lack 
of continuum emission indicates that such contributions 
are comparatively minor. 

Notably, however, the 24 /im emission (only) also 
shows the complete SNR shell of N49, including the west 




Fig. 10. — Spectral Energy Distributions (SEDS) for (a) N49 
and (b) N63A. Solid lines give the SED for the entire SNR; dotted 
lines give the SEDs for the three bright clumps in N63A. Filled 
squares give the 24 /an flux densities along portions of the SNR 
limb seen only at that wavelength; open squares with arrows give 
upper limits at other wavelengths. 



side which is largely undetectable in optical emission 
lines. The spectroscopy of N49 shows the presence of 
several lines within the bandwidth of the 24 /an image 
(AA w 4.7 /im), including an [Fe III] line at 23.0 /im, 
[Fe II] lines at 24.5 and 26.0 /im, and an [O IV] line at 
25.9 /im. The 25.9 and 26.0 /im lines are both consider- 
ably brighter than any of the lines which fall within the 
IRAC bands. Thus it is plausible that the source of the 
24 /im emission from the SNR shell may be primarily line 
emission. 

To quantify this estimate further, we summed the 
expected fluxes of the lines listed above, which, over 
the 11'.' 1 x 22'.'3 area of the slit, gave a total flux of 
2.2±0.2 x 10~ 12 erg citT 2 s" 1 . On the MIPS 24 /im im- 
age, we created an 11'.' 1 x 22'.'3 region at approximately 
the same sky location and orientation as the placement 
of the IRS Long-High slit. We summed over the sur- 
face brightness of this region of the MIPS image, with- 
out background subtraction (since the IRS flux estimates 
included both source and background), and obtained a 
total flux of 2.8±0.1 x 10" 12 erg cm" 2 s _1 . Thus, based 
on these rough estimates, line emission could be respon- 
sible for almost 80% of the MIPS 24 /im emission from 
N49. Note, however, that in addition to the statistical 
errors given, this estimate is subject to uncertainties in 
the calibrations of MIPS and IRS, particularly given the 
unknown background contribution to the IRS spectra. It 
is also quite possible that this estimate applies only to 
the comparatively IR-bright eastern region, and cannot 
be meaningfully applied to the rest of the limb, particu- 
larly the areas seen only at 24 /im. 

On the other hand, the flux density of thermal emission 
from hot dust is proportional to Q \B \(Tfi) , where B\ is 
the Planck function, Q\ is an efficiency factor for the 
infrared emissivity, and is the dust temperature. For 
metals and crystalline grains in the mid/far-IR, Q\ oc 
A~ 2 ; typical dust te mperatures from a young SNR are 
~ 100 - 130 K (e-g.- lStanimirovic et al"l t2005). Plotting 
this expected curve for the IRAC and MIPS wavebands, 
we see that the curve peaks in the 24 /im band. The curve 
predicts that the next highest flux density, in the 8 /im 
band, will be ~ 10 x fainter than that at 24 /im. Since the 
24 /im surface brightness of the northwest limb of N49 is 
relatively faint (~ 4 — 5<r above local background), it is 
quite plausible that emission ten times fainter would not 
be detected. 
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4.3.2. Origins of infrared emission in N63A 

For SNRs within H II regions, such as N63A, there may 
be sufficient far-ultraviolet radiation to destroy nearby 
PAHs. In fact, the shell (limb) of N63A is largely unde- 
tectable in any of the IRAC wavebands, which as we have 
noted above are most sensitive to PAH emission. At 24 
/im, however, the complete shell of N63A (also seen in the 
X-ray and radio regimes) is apparent. The radio and IR 
coincidence might result from a significant synchrotron 
component in the I R. However, extra polating the radio 
spectrum for N63A Ipickel et alJ lT9931 to 24 /im gives 
an expected flux density of 6 mJy at that wavelength, 
less than 1% of the observed flux density. Thus, we take 
the strong X-ray and IR correspondence as an indication 
that line emission again dominates in the shell at 24 /im. 
As with N49, this does not wholly preclude the presence 
of a hot dust component. However, the lack of optical 
emission again suggests a low-density environment (other 
than in the three-lobed nebula) and therefore that the re- 
gion lacks the necessary dense clumps for significant dust 
emission. 

In contrast, the three-lobed nebular region of N63A 
is apparent at all available wavelengths, most brightly in 
the western (photoionized) lobe. The spectral energy dis- 
tribution (SED) of the entire SNR, which is dominated 
by this central nebula, appears similar to that of N49. 
While it is plausible that these high-density regions in 
N63A may have similar emission properties to the high- 
density filaments of N49, more detailed spectroscopy is 
required to probe the dominant processes producing the 
infrared emission. When the lobes are considered sepa- 
rately, their SEDs are similar to one another overall; sur- 
prisingly, the SEDs of the southeastern (shocked) lobe 
and the western (photoionized) lobe have almost iden- 
tical spectral shapes, while the northeastern (shocked) 
lobe has an SED slightly divergent from those of the other 
two (Fig. 1 1 U|l . Given that the optical data suggest that 
the western lobe is an unshocked H II region while the 
eastern lobes are shocked by the SNR, it is surprising to 
find that the northeastern lobe, rather than the western 
lobe, shows the greatest difference from the other two 
lobes. 

4.3.3. Infrared and X-ray Comparisons 

One means of addressing the possible contributions of 
line versus continuium emission is to examine the rela- 
tive distribution of 24 /im emission c ompared to the X- 
ray em ission, as done, for example, bv lStanimirovic et al.l 
(2005) in the study of SNR IE 0102.2-7219. Following 
a similar procedure for N49 and N63A, we re-centered 
and re-gridded the Chandra X-ray image for each SNR 
to match the corresponding 24 /zm image. Each X-ray 
and 24 /im image was then normalized so that the peak 
emission associated with the SNR was set to 1.0. (For 
N49, the X-ray emission spatially corresponding to the 
location of SGR 0525-66 was not counted toward this 
"peak" emission.) The resulting normalization factors 
were 7 x 10~ 3 ct s^ 1 arcsecs~ 2 and 63 MJy sr" 1 for N49, 
and 5 x 10" 3 ct s _1 arcsecs -2 and 206 MJy sr" 1 for N63A. 
We then divided the 24 /xm images by the X-ray images 
in order to produce ratio maps (not shown), which give 
the distribution of 24 /im emission (with respect to the 24 
/im peak) compared to the X-ray emission (with respect 



to the X-ray peak). 

If we then consider an "infrared excess" to be present 
where these ratio maps show ratios greater than 1.0, and 
if we consider that regions dominated by thermal emis- 
sion from dust should show such an "infrared excess," 
then we can, by integrating over such regions, estimate a 
rough proportion of the total emission which, we pre- 
sume, is due to continuum emission from dust. For 
N49, roughly three-quarters of the flux could then be 
attributed to continuum emission! (Note that the IRS 
slit locations do fall roughly where the ratio maps show 
an excess of X-ray emission compared with IR, so this is 
not necessarily inconsistent with the overwhelming dom- 
inance of line emission in the IRS spectra.) For N63A, 
the "infrared excess" is largely seen on the western side 
of the SNR, including the three-lobed nebula; using the 
above criteria we would expect approximately half the 
emission to be continuum emission from heated dust. 

However, this som ewhat simplistic approach must be 
treated with caution. iStanimirovic et alTl)2005t) used the 
fact that that the strongest X-ray emission in IE 0102.2- 
7219 was associated with reverse-shocked ejecta in that 
young SNR, a region which would be expected to gen- 
erate strong line emission in many wavelength regimes. 
In contrast, N49 and N63A are both older remnants, in 
which the ejecta is probably intermixed with swept-up 
ISM to a large extent. For both these SNRs, some re- 
gions identified as having an "infrared excess" are also 
sites of significant optical line emission - including parts 
of the southeastern limb of N49 and the "three-lobed 
nebula" of N63A. IR emi ssion associated with optically- 
emitting was suggested bv lArendt et al.l (^992) to include 
significant infrared line emission! Other regions, such 
as the western side of N49 and the southwestern limb 
of N63A (the "X-ray loops"), probably lack significant 
optical emission because that region of the SNR is ex- 
panding into a low-density environment. This raises the 
point that we would not expect the high-density clouds 
required to produce substantial thermal dust emission in 
such regions. 

4.3.4. IRAC colors 

iReach et all (|2006l) . in their study of Galactic SNRs, 
were able to measure fluxes for a number of those SNRs 
in all four IRAC bands. Based on measurements and sim- 
ple models, the authors developed an IRAC color-color 
diagram showing the regions of "color space" expected 
for different emission mechanisms, including PAH, ionic, 
molecular and synchrotron emission. They found lit- 
tle evidence of synchrotron emission for their detected 
SNRs; and state that "observations of infrared emis- 
sion from SNRs to date have shown little or no evidence 
of significant dust emission (PAH or continuum) within 
the wavelengths of the IRAC bands." From their sam- 
ple, they identified regions within the SNRs with colors 
typical of emission from molecular lines, ionic lines, and 
PAHs. A number of their detected SNRs showed signifi- 
cant molecular-line emission; several of these have other 
signs of molecular cloud interaction, such as maser emis- 
sion. 

N49 and N63A are the only SNRs in this sample 
for which all four IRAC bands are available. Using 
N49, we determine colors of I(3.6/im)/I(5.8/im)=0.25 
and I(4.5/xm)/I(8/W)=0.23. These colors place N49 be- 
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TABLE 5 

SNR SIZES AND AGES 



SNR 


Angul 


ar 


Linear 


Estimated 


Age 


Name 


Radii 


a 


Radii 


Age 


reference 


NHL 


35" x 


30" 


8.8 x 7.5 pc 


7000-15,000 yr 


Williams et al. (1999) 


N44-SNR 


134" x 


122" 


33.5 x 30.5 


~18,000 yr 


Chu et al. (19931 


N49 


40" x 


40" 


10.0 x 10.0 pc 


~6600 yr 


Park et al. (2003) 


N63A 


37" x 


37" 


9.3 x 9.3 pc 


2000-5000 yr 


Hughes et al. (1998) 


N206-SNR 


95" x 


95" 


23.8 x 23.8 pc 


23,000-27,000 yr 


Williams et al. 72005) 


N157B 


68" x 


58" 


17.0 x 14.5 pc 


~5000 yr 


Wang fc Gotthelf (1998) 



a Size is based on the semimajor and scmiminor axis of an ellipse fit to the 
optical extent of the SNR; except for N49 and N63A, where the X-ray extent 
is used because some of the SNR is not seen in optical. 



tween the "molecular shocks" a nd "ionic shocks" sec tions 
of the color-color diagram of Rea ch et alJ l)2006t Fig. 
2), somewhat closer to the area where molecular shocks 
predominate. These colors are thought to represent "a 
mixture of molecular and ionic shocks," and are seen 
in such remnants as, e.g., CTB 37A and G348.5— 0.0. 
We see from the ionic and molecular lines identified by 
Spitzer IRS spectroscopy (Table0J t hat this is indeed an 
accurate characterization for N49. iReach et alJ l)2006f) 
notes that both CTB 37A and G348. 5-0.0 are associ- 
ated with OH 1720 masers, which strongly suggests that 
these remnants are interacting with molecular clouds, as 
N49 is thought to be. All lobes of the three-lobed nebula 
in N63A have colors (Table [3Jl that place them within 
the "molecular shock" portion of the diagram, suggest- 
ing that shocked molecular gas dominates at the near-IR 
wavelengths. 

5. CONCLUSIONS 

The point that line emission is often a substantial 
contributor to the infrared emi ssion from SNRs has 
been raised for so me time (e.g . . iGraham et alJ fl987bt 
lArendt et all IT992T: loliva et alJ 119971 Il999aj) . For in- 
stance. lArendt"et al.l l|1992f) found that infrared line 
emission from shocked gas contributes "a significant frac- 
tion" to one component of the infrared emission - the 
component generally traced by the presence of optical 
emission. (Note that dust emission was also thought 
to contribute to this component, and that a second IR 
component, traced by an excess of IR to X-ray emis- 
sion, was inferred to be gene rated only by co llisionally 
heated dust.) More recently, IStanimirovic et all l)2005fl 
find that for the Small Magellanic Cloud SNR IE 0102.2- 
7219, emission from the [O IV] line may contribute up to 
60% of the 24 /im emission. 

Our observations of these remnants have primarily 
shown infrared emission, where present, arising from 
optically-bright areas of the SNRs. The implication is 
that line emission is a significant contributor to the in- 
frared emission of these re mnants. The IRAC-band ra- 
tios of IReach et all l|2006|) for the optical/IR-bright ar- 
eas in N63 and N49 also indicate substantial contribu- 
tions from ionic and molecular lines. This suggestion 
is greatly strengthened by the infrared spectroscopy of 
the optically-bright limb of N49, where the line contri- 
butions appear to overwhelmingly dominate the infrared 
flux density. 



Perhaps the most interesting (and mysterious) result is 
that for our two clearest detections, N49 and N63A, the 
outer limb of the SNR (in both cases, extremely faint in 
optical emission lines) is dete cted clearly at 24 um but 
not in the other wavebands. IStanimirovic et alJ l)2005f) 
find a similar result for IE 0102.2-7219; the SNR is de- 
tected clearly at 2 4 /im , but not at 8 or 70 /im. Similarly, 
iBorkowski et all £2006), in a study of Type la SNRs in 
the Large Magellanic Cloud, find that none of their four 
SNRs are detected in the IRAC bands; two are detected 
at 24 and 70 /im, and two at 24 /im alone. 

The comparative brightness of emission from these 
SNRs in the 24 /im band, compared to the other IR 
bands, may be due in part to the line emission from 
[O IV] and [Fe II], particularly in regions of bright op- 
tical emission. However, we are far from able to rule out 
a secondary dust component even in these regions. The 
areas of strong 24 /im emission, weak optical emission, 
and moderate X-ray emission may point to a very dif- 
ferent physical case, where IR emission from collisionally 
heated dust may yet predominate. Our future spectral 
mapping of N49 and N63A with Spitzer's IRS are ex- 
pected to elucidate this situation. 

The potentially significant contributions of infrared 
line emission from at least some SNRs has serious im- 
plications for the calculations of the dust content within 
those SNRs, and the use of such calculations to constrain 
models of dust production and destruction. For example, 
dust mass estimates based on the 24 /im Spitzer flux or 
the 25 /im IRAS flux may overestimate the amount of 
dust due to overestimation of the flux resulting from in- 
clusion of line emission. For some remnants, e.g. where 
nonradiative shocks dominate, this difference may be 
only slight, but for highly radiative SNRs, particularly 
those interacting with dense ambient material, the effect 
may be significant. Likewise, fits to SEDs may be un- 
physical if line contributions in one or more bands are 
not properly taken into account. 

This work is based on observations made with the 
Spitzer Space Telescope, which is operated by the Jet 
Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA. Support for this work 
was provided by NASA through JPL/Caltech award 
JPL-1264494. RMW is supported by NASA under grant 
NNG05GC97G through the LTSA program. 
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